Observational constraints on unified dark matter including Hubble parameter data 
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We constrain a unified dark matter (UDM) model from the latest observational data. This model 
assumes that the dark sector is degenerate. Dark energy and dark matter are the same component. 
It can be described by an affine equation of state Px = po + apx- Our data set contains the 
| newly revised H(z) data, type la supernovae (SNe la) from Union2 set, baryonic acoustic oscillation 

(BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 (DR7) 
galaxy sample, as well as the cosmic microwave background (CMB) observation from the 7-year 
Wilkinson Microwave Anisotropy Probe (WMAP7) results. By using the Markov Chain Monte 
Carlo (MCMC) method, we obtain the results in a flat universe: r2A=0.719lo;o3os(lo-)lo.o458( 2cr )> 
a=1.72t^(lcr)t^(2o-)(xl0- 3 ), fi fl / i 2 =0.0226+J!;°°^ (la)lg;°°^(2cr). Moreover, when consid- 
ering a non-flat universe, n A =0.722l°;°^(lcr)±°;°*^(2cr) ) a=0.242l°;^(lo-)l^°(2cr)(x 10" 2 ), 
f7 6 /i 2 =0.0227i°°°^(la)l|;'°°^(2cr), n k =-0.l94+f° a 2 5 {lo)t\ :£r(2<r)(x 1(T 2 ). These give a more 
stringent results than before. We also give the results from other combinations of these data for 
comparison. The observational Hubble parameter data can give a more stringent constraint than 
O f SNe la. From the constraint results, we can see the parameters a and Qk are very close to zero, 

which means a flat universe is strongly supported and the speed of sound of the dark sector seems 

5h " 

■ to be zero. 

c/2 
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Cosmic acceleration is one of the most striking discoveries in modern cosmology [1] . Many works have been done to 



INTRODUCTION 



explain this unexpected phenomenon. There are different mechanisms for the universe acceleration. The most popular 
idea is assuming a new component with negative pressure known as dark energy. Other mechanisms contain various 
modifications of gravitation theory Q and violation of cosmological principle Q. There are all kinds of dark energy 
models including quintessence Q], holographic dark energy quintom phantom generalized Chaplygin gas 
model [8] and so on. Among these, the simplest and most successful model is ACDM model Q while it has to face 
some theoretical problems. For example, the "coincidence" problem and the "fine-tuning" problem [loll. From this 
standard model, the flat universe consists of ~ 4% baryons, ~ 23% dark matter, ~ 73% dark energy However, 
we still know little about dark matter and dark energy jl2j. All we can feel about them are their gravitational effects. 
Observationally, they arc degenerate. Based on these, a unified dark matter model was proposed which assumes that 
dark matter and dark energy are from a same dark component Since the equation of state of the unified dark 
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component is unknown, we assume that it has a constant speed of sound c 2 s — a which is equivalent to the affine form 
[l^j . In order to break the degeneracies, different kinds of observational data are used to constrain cosmological model, 
among which SNe la, CMB, and BAO are usually used to constrain cosmological parameters 
observations are included. For example, the gamma-ray burst [16| , cluster gas mass fraction [17 



15| . Rece ntly , other 
, lensing [18( and so 



on. But these processes need to integrate the Hubble parameter to obtain the distance scale. The more information 
about H(z) can't be embodied. Therefore, it is meaningful to study the H(z) data directly. The expression of the 
Hubble parameter can be written in this form 

1 + z at 

which depends on the differential age as a function of redshift. Based on Jimenez et al. [2^], Simon et al. 21 1 used the 
age of evolving galaxies and got nine H{z) data. Many works based on these nine data was done 22J. Recently, Stern 



et al. 



et al. 25 1 revised these data at 11 redshifts from the differential ages of red-envelope galaxies. Moreover, Gaztahaga 



26[ took the BAO scale as a standard ruler in the radial direction, obtained three more additional data. 



In this paper, we use the Markov chain Monte Carlo (MCMC) method to give constraints on the unified dark 
matter model from the latest data including the observational Hubble parameter data. In section 2, we give a brief 
introduction of the UDM model. In section 3, we introduce the observational data set we use. The full parameter 
space using different combinations of data is exhibited in section 4. Finally, we give a conclusion in section 5. 



2. UNIFIED DARK MATTER WITH CONSTANT SPEED OF SOUND 

In this section, we give a brief introduction to unified dark matter model. For details of this model, please refer to 

Considering the universe is described by the Friedmann-Robertson- Walker metric and assuming a single UDM 
besides baryons and radiation, the Friedmann equation is expressed as 

H(a) = H E(a), (2) 

where Hq represents the Hubble constant today, and 

E{af = n r (a) + Sl b (a) + Sl x (a) + fi fc (o), (3) 

where f2j(a) = pi(a)/p c , p c is the critical density today. The subscripts r, b, k and X represent for radiation, baryons, 
curvature and UDM component respectively. Each component satisfies the energy-momentum conservation equation 

Pi = -W{p i + P i ). (4) 

For the radiation and baryons, the equations of state are 1/3 and 0, respectively. The conservation equation gives 
Qb(a) oc a~ 3 and Q r (a) oc a~ 4 . 

For the UDM, it is considered as an approximation to an barotropic fluid. Since the equation of state Px = Px(px) 
is unknown, we assume a constant speed of sound 



dPx/dpx - «, 



(5) 
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which leads to the 2-parameter affine form 



Px —Po + apx- 



(6) 



Combined with the conservation equation, we can regard this component as an effective cosmological constant plus 
a generalized dark matter which has a constant equation of state a. The density evolves as 



where O m = flxo ~ stands for the generalized dark matter, and stands for the effective cosmological constant 
component, p\ = 0. 

The Friedmann equation now can be expressed as 



where fl r , fib and Qfc represent dimensionless radiation, baryon matter and curvature today. 51a and a are the 
parameters we introduced. We notice that when a = 0, the UDM model are equivalent to ACDM. In 
UDM model, we can equivalently think that the generalized dark matter may has a non-zero equation 
of state. Naturally, because the property of the effective cosmological constant component is the same 
as the one in ACDM, when the property of generalized dark matter is equivalent to the usual dark 
matter in ACDM, UDM model totally becomes ACDM. 

The speed of sound characterizes the perturbations to the dark matter density and pressure which 
are related to the microphysics. Since we know litter about dark sector in our universe, the equation 
of state is introduced to describe the property of dark sector. However, it can not give us the 
information about the microphysics and other properties. The speed of sound can give information 
about the internal degrees of freedom. For example, quintessence has c s = 1. If we can confirm the 
sound speed that is below the speed of light, we can get further degrees beyond a canonical, minimally 
coupled scalar field. If the sound speed is small, the perturbations can be detected on a proper scale. 
This perturbations will change the gravitational potential, thus affecting the geodesic of photons. As a 
result, the cosmic microwave background will be changed through the Integrated Sachs- Wolfe (ISW) 
effect. In order to show the differences between UDM mdoel and ACDM which corresponds to a = 0, 
we give the CMB temperature angular power spectrum in Fig. 7. Our calculations are based on 



-3(l+a) 



(7) 



E(a) 2 = ft r a- 4 + n b a- 3 + fl A + (f - fl b - O r - rj A )a _3(1+a) + n k a~ 2 , 



CAMB [19]. 



3. 



CURRENT OBSERVATIONAL DATA 



Now, we introduce the data set using to constrain the UDM model. 



3.1. Observational Hubble parameter data 



The observational Hubble parameter as a function of rcdshift has been used to constrain cosmological parameters 
23 1 . For the observational data, we choose 11 data from [25 1. These data were obtained through the measurement of 
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differential ages of red-envelope galaxies known as " differential age method" . The aging of stars can be regarded as 
an indicator of the aging of the universe. The spectra of stars can be converted to the information of their ages, as 
the evolutions of stars are well-known. Since we can not observe stars one by one at cosmological scales, we take the 
spectra of galaxies which contain relatively uniform star population. The surveys including the Gemini Deep Deep 
Survey 24 1, VIMOS-VLT Deep Survey and SDSS Q are related to obtaining data. We consider the red-envelope 
galaxies which are massive with fairly homogeneous star population. Moreover, the data can be obtained from the 
BAO scale as a standard ruler in the radial direction known as "Peak Method" [26j. Though this method, we have 



three additional data [26(: H(z = 0.24) = 79.69 ± 2.32, H(z = 0.34) = 83.8 ± 2.96, and H(z = 0.43) = 86.45 ± 3.27, 
which are independent. Plus the measurement of Hubble constant we totally have 15 data now. The x 2 value of 
the H(z) data can be expressed as 

2 ^ [H{ Zi ) - H obs (z t )} 2 

■;=i 1 

where <ii is the la uncertainty of the observational H(z) data. 



3.2. Type la supernovae 

Since the SNe la fist revealed the acceleration of the universe, it has been a well-established method for probing 
dark energy. 



Recently, the Supernova Cosmology Project (SCP) released the Union2 data set including 557 samples [28]. It 
covers a redshift region 0-1.4. These data have been widely used to investigate dark energy models 29]. The distance 
modules is given by 

^theory(z) = 5 log lf) [d L (z)} + fl 0l (10) 

where cIl(z) is the luminosity distance. It is defined as 

r , r z Ay' ' 

vW/ 4h ■ (ii) 

The marginalized nuisance parameter [30( for x 2 is expressed as 



d L (z) = (1 + z)r(z), r(z) = -— ==sinn 

#ov \\lk\ 



E(z>) 



where A = £f 7 (/x data - fi thcor y ') 2 / 'of , B = ^f 7 (^ data ~ ^ hcmy )/c 2 , C = £)f 7 1/of , a % is the la uncertainty of 
the observational data. 
Note that the expression 

xh c = A-^, (13) 
which coincides to Eq. (|12[) up to a constant [31], is used in this paper. The constraint results are the same. 



5 



3.3. Cosmic microwave background 



For CMB, the shift parameters R and l a are important indicators which are related to the positions of the CMB 
acoustic peaks. These quantities are affected by the cosmic expansion history from the decoupling epoch to today. 
The shift parameter R is given by 32 1 

'' ' dz 



1/2 



(14) 



/o E{z)_ 

where S7 m o = Q m (hereafter), since we have defined the generalized dark matter. The redshift corresponding to 
the decoupling epoch of photons is given by [331 ] 



= 1048[1 + 0.00124(fM 2 )-°- 738 (l + gi(n m0 h 2 ) 92 )}, 



(15) 



with 5 i = 0.0783(fVi 2 )~ - 238 (l±39.5(lM 2 )~ a763 r\ 52 = 0.560(1 + 21.1(n 6 ft 2 ) 1 - 81 )- 1 . Moreover, the acoustic scale 
is related to the distance ratio. It is give by 



la = 7T 



where the comoving sound horizon at photo-decoupling epoch is expressed as 



r s (z*) = i? 1 / c s (z)/E(z)dz. 



From WMAP 7 measurement, the best- fit values of the data set is given by [3f 









( 302.09 ±0.76 \ 


PCMB = 


R 




1.725 ±0.018 




J 




^ 1091.3 ±0.91 j 



where APcmb = Pcmb ~ Pcmb, and the corresponding inverse covariance matrix is 

/ 2.305 29.698 -1.333 \ 



'CMB 



29.698 6825.270 -113.180 
\ -1.333 -113.180 3.414 J 



With these, we can calculate the x 2 by 



35] 



Xcmb — APcmbCcmb 1 APcmb- 



(16) 



(17) 



(18) 



(19) 



(20) 



3.4. Baryon acoustic oscillation 



The Baryon Acoustic Oscillations we use measure the distance-redshift relation at two redshifts. One is in the 
clustering of the combined 2dFGRS and SDSS main galaxy samples at z = 0.2, the other is in the clustering of the 
SDSS luminous red galaxies at z = 0.35. Calculated from these samples, the observed scale of the BAO are analyzed 
to constrain the form of the distance scale. The distance scale of BAO is defined as 

1/3 



D v (z) 



n k H(z) 



sum 



dz' 
o HW) 



(21) 
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The peak positions of the BAO depend on the ratio of D v (z) to the sound horizon size at the drag epoch. The 
redshift Zd corresponding to the drag epoch when baryons were released from photos is given by 

1291(f2 m0 /l 2 ) - 251 r , , »,2\&2Y1 fool 

Zd = [l + 0.659(n m0 ^)°-^] [(1 + W» h ) ')]» ( 22 ) 



where h = 0.313(^ m0 /i 2 )-°- 419 [l + 0.607(fi m0 /i 2 ) a674 ]" 1 and b 2 = 0.238(fi m0 /j 2 ) - 223 [36j . In this paper, we choose 
the measurement of the distance radio (d z ) at z = 0.2 and z — 0.35 37:]- This quantity is given by 

* - wiry (23) 

where r s (zd) is the comoving sound horizon. From the SDSS data release 7 (DR7) galaxy sample, the best-fit values 
of the data set (do.2, do. 35) [37j 



p 

A matrix 


' do.2 1 


-l 




1 do.35 j 





0.1905 ±0.0061 
0.1097 ±0.0036 



(24) 



where the corresponding inverse covariance matrix is 



_ x I 30124 -17227 \ 

Cmatrix = • (25) 

I -17227 86977 / 
The x 2 value of this BAO observation from SDSS DR7 can be calculated as 

XBAO = APmatrixCmatrix 1 AP matr i x . (26) 

where AP ma t r j x Pmatrix Pmatrix- 



4. CONSTRAINT RESULTS 



The statistics we use to constrain parameters is the usual maximum likelihood method of \ 2 fitting with the Markov 
Chain Monte Carlo (MCMC) method. Our code is based on CosmoMCMC [38]. The corresponding convergence of 
the chains R-l is set to be less than 0.003 which can guarantee the accuracy well. We combine H(z), SNe la, CMB 
and BAO data to give a global fitting on determining the parameters of UDM model. The total \ 2 1S given by 

Xtot =X 2 H+ XSNc + XCMB + XBAO- ( 27 ) 

Fist, we consider a flat universe. In order to embody the H{z) data, we choose three combinations of data 
sets and show the 1-D probability distribution of each parameter (Tl^h 2 , Cl\, a) and 2-D plots for parameters 
between each other. Fitting results from the joint data of i?(z)+CMB+BAO are given in Fig 1, the best- 
value for a is a=1.6lt4;5g(lcr)lgg2(2cr)(xl0~ 3 ). For comparison, we give the results from SNe±CMB+BAO in 
Fig 2. The best-fit value is a=1.81+g '^(ler) j^'g^er^x 10~ 3 ). We can see the observational Hubble parameter 
data gives a more stringent constraint than SNe la in this model. In Fig 3, the total combinations of them re- 
sult in a= lJ2tl 9 7 2 9 (l<7)t 5 7 il(2a)(xl0- 3 ). From the result, the parameter a which stands for the speed of sound 
of the unified dark sector seems to be vanishing or slightly larger than 0. So ACDM model is strongly sup- 
ported. We further consider a non-flat universe with curvature. For i/(z)+CMB-(-BAO in Fig 4, the results 
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FIG. 1: The 2-D regions and 1-D marginalized distribution with the 1-a and 2-er contours of parameters Qbh 2 , a, Qa in a flat 
universe, for the data sets (z)+CMB+BAO. 

are a=2.2lt 7 7 '^(la-)t\o^(2a)(xlO- 3 ), il k =-0.208tf^(la)tl\^(2a-)(xl0- 2 ). Using SNc+CMB+BAO in Fig 5, 
a=0.662^ 2 8 ° 2 (lcr)t^^(2cr)(xl0" 3 ), fi fc =-1.4lt|;^(lcr)t^^(2cr)(xl0" 2 ). This results are consist with the situa- 
tion in a flat universe. The H(z) data also gives a much more stringent constraint on Slfc. The total joint data in 
Fig 6 give a=2.42^;^(lcr)+ 1 J;°°(2cr)(xlO" 3 ), n k =-0A94tlf 5 (l<y)tli 7 '(2a)(xl0- 2 ), which suggest both the speed 
of sound and the curvature are very little. Our results are much more stringent than previous work [l4j |. Our results 
seem to be consistent with a = using the data sets we select. Does that mean there is no need to 
introduce a new parameter? We think although the UDM model is contrived, the sound speed of the 
unified dark matter is related to the microphysics of the dark component while equation of state does 
not. The sound speed leads to the inhomogeneities of dark matter, thus affecting CMB and matter 
power spectra. With the development of technique, more information in CMB, matter power spectra 
and so on may distinguish UDM model from ACDM. Moreover, we assume a constant sound speed, 
the whole effective value is very small, however, if when we consider the evolving speed of sound, the 
results may be changed. 



5. CONCLUSIONS 



In this paper, we constrain the unified dark matter with constant speed of sound from the latest data: the newly 
revised 15 Hubble data, type la supernovae (SNe la) from Union2 set, baryonic acoustic oscillation (BAO) ob- 
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FIG. 2: The 2-D regions and 1-D marginalized distribution with the 1-a and 2-a contours of parameters Qbh 2 , a, Oa in a flat 
universe, for the data sets SNe+CMB+BAO. 





Flat Universe 




#(z)+CMB+BAO 


SNe+BAO+CMB 


#(z)+SNe+CMB+BAO 


Oa 
a 


0.717lH^(l-)tH^(2-) 
1.6ll|;f 8 (l C T)i^(2a)(xlO- 3 ) 

0.02261™^ (l^S^) 


0.719^-° 2 3 f 3 (la)« : »« 2 2 (2a) 

l.8ii 8 8 ;f (ia)^i; 9 6 9 9 (2<T)(xlo- 3 ) 

0.0227l°;^(la)lg : ^(2a) 


0.719lS:SS(la)lS:^(2a) 
1.72+ 3 4 .? 2 (la)i^(2a)(xl0- 3 ) 
0.0226t»;S-(la)lg : »; 6 5 (2a) 


Xmin 


11.592 


532.142 


542.354 



TABLE I: The best-fit values of parameters Oj,/i 2 , a, Qa for UDM model in a flat universe with the 1-a and 2-a uncertainties, 
as well as Xmin, for the data sets J7(z)+CMB+BAO, SNe+BAO+CMB, and #(z)+SNe+BAO+CMB, respectively. 



servation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 (DR7) galaxy sample, as well 
as the cosmic microwave background (CMB) observation from the 7-year Wilkinson Microwave Anisotropy Probe 
(WMAP7) results. We consider two case. In a flat universe, the constraint results are flA=0.719too305(lc0^g'o45g(2er), 
a=1.72±l%l(la)t 5 7 il(2o-)(xW- 3 ), n b h 2 =0m26toiou(^)to^o\U 2a )- In a non-fiat universe, the constraint re- 
sults are O A =0.722l°;^f 7 (la)^;^ 9 4 (2 ( T), a=0.242^(lG-)t]:^(2o-)(xlO^), n b h 2 =0.0227t Zl(^)^l(2cT), 
flk=— 0.194t 2 /g 2 l (lcr)i 2 , 57(2<7)(xl0~ 2 ). This constraint results are much stringent than previous work 14|, which 
suggest both the speed of sound of the UDM and the curvature are very little and a flat ACDM model is still 
a good choice. In order to know the effects of the observational Hubble data, we also consider the combinations 
i?(z)+CMB+BAO and SNe+BAO+CMB. From the results, the H(z) data gives a better constraint especially for a, 
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FIG. 3: The 2-D regions and 1-D marginalized distribution with the 1-a and 2-cr contours of parameters Qbh 2 , a, Oa in a flat 
universe, for the data sets i/(z)+SNe+CMB+BAO. 





Non-flat Universe 




#(z)+CMB+BAO 


SNe+BAO+CMB 


J ff(z)+SNe+CMB+BAO 


a 
Q b h 2 


0.721±g : gg*(la)lg : »«(2a) 
2.2ll^ (la)+^(2a)(xl0- 3 ) 

-0.208t?;^(la)^;^(2a)(xl0- 2 ) 


0.734i° ;^(la)l :?S 9 (2 f r) 
0.662^ 2 8 ° 2 (la)+f 7 6 3 8 3 (2a)(xl0- 3 ) 

0.0220l°;^(la)l°;S(2a) 
-1.4lt 6 4 | 5 3 (la)l 9 6 ;^(2a)(xl0- 2 ) 


o.722±™(i<7)i°o: 63r(2*) 

2.42^(la)^J;»° 9 (2a)(xl0- 3 ) 

0.0227l°;r^(l-)lH^8(2-) 
-0.194l 2 ;^(la)i 2 2 .I 5 7 (2a)(xl0- 2 ) 


Xmin 


11.578 


531.870 


542.309 



TABLE II: The best-fit values of parameters Qbh , a, Qa, fife for UDM model in a non-flat universe with the 1-a and 
2-cr uncertainties, as well as xLin, for tne data sets (z)+CMB+BAO, SNe+BAO+CMB, and #(z)+SNe+BAO+CMB, 
respectively. 

fife and fib/i 2 than SNe la in this model. We think the H(z) data will paly an important role in cosmology in the 
future Q. 
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FIG. 4: The 2-D regions and 1-D marginalized distribution with the l-a and 2-a contours of parameters Q.bh 2 , a, Oa, fi* in a 
non-flat universe, for the data sets H (z)+CMB+BAO. 
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